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Abstract 

The profile functions of the SU(3) Skyrme soliton are investigated for the octet, decuplet, and 
antidecuplet baryons by the mean field approach. In this approach, the profile functions are affected 
by the spatial rotation, the flavor rotation, and the flavor symmetry breaking. The solitons are 
stable only in the restricted areas of the parameter space for each multiplet. When the flavor 
symmetry breaking is large, the area for the antidecuplet is narrow compared to those for the octet 
and decuplet. The parameters are determined by the baryon mass spectrum, and the deformation 
of the soliton has sizable effects on the masses. 
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I. INTRODUCTION 



Diakonov, Petrov, and Polyakov [jj made a detailed prediction for the masses and the 
decay widths of the antidecuplet (10) baryons in the framework of the Skyrme sohton 
(Skyrmion) model . Following their work, an experimental discovery of the light- 

est state of 10, namely 0^(1540), was reported by the LEPS collaboration Is]. 6"^ has 
strangeness = +1 and should contain at least one s quark. It is called an exotic baryon 
or a pentaquark, because the minimal number of the quarks is five from the charge and the 
strangeness. Although later many experiments confirmed this finding, several experiments 
did not observe Lists of these published experiments and detailed discussion of their 
results are presented in Refs. g], 7|. 

the diquark mod- 

the QCD sum rules [iSj], and the lattice QCD [14]. These works are reviewed in Ref. 
We are interested in the descriptions of the 10 baryons by the soliton 



Theoretically, there are many works based on the Skyrme model 
els qI, the chiral bag model 10], the MIT bag model 11|, the constituent quark mode^ 
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3, 
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181] in the 



231 ]. Now, there are two major approaches to quantize 



SU(3) Skyrme model [19| 
the soliton. First is the Callan-Klebanov approach 2J], in which baryons appear as kaon- 
SU(2) Skyrmion bound states, and the isospin rotation of the soliton and the fluctuations 
of the kaon field are quantized. The bound states change according to the baryon states. In 
particular the Wess-Zumino term acts as a repulsive force on the S" = +1 states; its strength 
is strong enough to remove all bound states {3] and all resonances 26] for the standard 
values of the parameters. However, recently, Itzhaki et al. [lol applied this approach to 
the exotic baryons and found the kaon bound states of S* = +1 by using a large kaon 
mass (~ 1 GeV). ^ 

Second is the rigid rotator approach (RRA) [3], in which the shape of the soliton is 
common to all baryon states and the rotation of the soliton in flavor space is quantized. Then, 
the baryons emerge as the rotational states of a rigid soliton. From early papers 2l|, |22] 
on the SU(3) Skyrme model, it was pointed out that this approach reproduces not only the 
octet (8) and decuplet (10) baryons but also the antidecuplet (10) baryons as the low lying 
spectrum. The 10 baryons have the spin and the parity = 1/2"*" in this a ppr oac 



n. 



However, a limit of the applicability of RRA has been pointed out 
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29j in the 



SU(2) Skyrme model. The shape of the soliton changes because of the centrifugal force of 
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the rotation, and the large spin of the baryon leads to the instability of the soliton due to 
the spontaneous emission of the real pion from the soliton. In the SU(3) Skyrme model, the 
low lying multiplets (8, 10, and 10) seem free from the limit of the applicability of RRA l| 
due to their small spins. However, the situation depends on the baryon states. The rotation 
emerges in the strangeness direction simultaneously and pushes the shape of the soliton out 
further. In addition, if the shape of the soliton is affected by the strangeness degrees of 
freedom, it would shrink because of the large meson mass. Therefore, there is a possibility 
too that the flavor symmetry breaking cancels the deformation caused by the rotation. We 
consider that this possibility should be investigated particularly. 

Furthermore, Itzhaki et al. 1^ and Cohen |l7| pointed out that in a large number of the 
color {N^ expansion, the mass differences between the 10 and 8 baryons scale as A*"^. This 
means that RRA for multiplet 10 is not consistent with the large N^. expansion. Since the 
above mentioned deformations of the classical soliton are formally sub-leading effects in the 
expansion, we are interested in whether the effects could be practically negligible in 10. 

In this paper, we formulate a mean field approach to include the effects of the rotation 
and the symmetry breaking into the shape of the soliton. In this meaning, we modify RRA. 
In addition, we study numerically the soliton solutions derived from our approach, and find 
the input parameters that keep the soliton stable and reproduce the baryon mass spectrum. 

In Sec. ini the SU(3) Skyrme model and its collective quantization are reviewed, and the 
mean field approach for the soliton is introduced. In Sec. IIIIl the stability conditions of the 
soliton solution are explained, and the numerical solutions are displayed. In Sec. HVl the 
input parameters and the resultant baryon mass spectrum are given. Finally, in Sec. |V] we 
summarize the results. 
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II. SU(3) SKYRME MODEL 



A. Model and the collective coordinate quantization 



The effective action [30|] we take here is given by 



r — r5 + TsB + ^wz, 



(1) 




(2) 



(3) 



(4) 



where is the pion decay constant, e is the Skyrme parameter, U{x) is the SU(3) unitary 
matrix representing the pseudoscalar mesons (tt, K, rj), mT,^K,r} are the masses of (vr, K, 
7]), Ag is the 8th component of the Gell-Mann matrices (/U = 1,2, .. . ,8), and is the 
number of color degrees of freedom. The Wess-Zumino term Twz is given as an integral over 
the five-dimensional disk D^, the boundary of which is the compactified space-time 5*4. The 
symmetry breaking mass term Tsb contains only two masses rri-^^K because of the quadratic 
sum rule: + 3m^ — Amj^ = 0. In this paper, we choose (e, /t^, 111^,, itik) as adjustable 
parameters. 

The effective action ([T]) admits a classical static soliton solution under the hedgehog 
ansatz embedded in the SU(2) subgroup: 



where r = |r|, = Xj/r, and F{r) is the profile function of the soliton. The baryon number 
one solution is subjected to the boundary conditions 





F(0) = TT, F(oo) = 0. 



(6) 



( — I 



We postulate the cranking form [3] of the time dependent meson field: 



U{x) = Ait)Uoir)A\t), 



(7) 
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where A{t) describes the adiabatic col' 



Using the standard method 



19 



20 



21 



ective rotation of the system in SU(3) flavor space. 
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23| to quantize the motion on the SU(3) group 



manifold, we obtain a dimensionless quantized collective Hamiltonian H and a first class 
constraint on the 8th generator Rg of SUr(3): 



H 



Rft 



+^C2(SUR(3)) + |[l-Dg)(A) 



2v^' 



(8) 
(9) 



where C2(SUr,(3)) and C2(SUr,(2)) are the Casimir operators of SUr(3) and SUr(2) respec- 
tively. In addition, 

(10) 



DfJ{A) = -TT{XsA^XsA) 



and 
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— J dpp^ sin^ F f F'2 + 1 + — ^ 



4:71 / dpp^ sin' 



F fF 
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4 2,= 

dpp^{l — cosF), 



(11) 
(12) 
(13) 
(14) 



where p = ef-„r, rfiT^^K = i^tt,k /{ef-K), and F' = The Hamiltonian H explicitly depends 
only on e and friT^^K- The Hamiltonian, the classical soliton mass, and the symmetry breaking 
with the physical unit (MeV) are given by H = ^H, Mq = ^M^, and 7 = ^7 respectively. 
The moments of inertia with the physical unit (1/MeV) are given by = a'^/le^fj^) and 

The state function of the baryon B is labeled as 



B 




(15) 



where (J, J3), (/, I3), Y, and are the eigenvalues of the spin, the isospin, the hypercharge, 
and the right hypercharge, respectively. The right hypercharge is given by Y/^ = 1 due to 



constraint ([9]). Equation (fT5l) is obtained by solving the following eigenvalue equation, 



H'^B = Eb^b, 



(16) 



where Eb is the dimensionless energy eigenvalue of the baryon state '^b- To solve this 
equation, we use the Yabu and Ando method {sil, in which Eb is given by 



E 



B 



Quantity £sb is the dimensionless eigenvalue of 



C2(SUR(3)) + ^(l-Dg)(A)) 



B — £-SB^ B- 



(17) 



(18) 



B. Mean field approach to the baryon states 

To solve the eigenvalue Eq. f|T6|) and obtain the baryon states j^, 4, 
about the profile function F{p) in Eqs. (fTT!) - (fT4l) . Here, we define an equation of motion for 

F{p) as 



3l|, one should know 



B 



6F{p) 



0, 



(19) 



where T-Cb is a classical Hamiltonian for each baryon B. In RRA j^, Hb = Mo? and F is 
not affected by the rotation and the symmetry breaking at all. Therefore, F is common to 
all baryons. 

In this paper, we adopt the following mean field Hamiltonian 



32 



33 



34, 



nB = {^^B''\H\K') 



(20) 



where is an eigenstate of H without the SU(3) symmetry breaking (7 = 0), and the 

state is represented by the SU(3) D function corresponding to the baryon B. 

Our classical Hamiltonian Hb includes the influence of the flavor symmetry breaking of 
the first order in powers of {ml^ — m^) and the influence of the rotation of order l/Nc'm the 
large expansion. The specific expression of Tis is given by 



Ti^ij = Mo + - [a'^' + r j + -gB7, 



(21) 
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TABLE I: SU(3) representation {p,q), spin J, J(J+1) in 


Eq. ([22]), andC2(SUR(3))- 


^(^+1)-! 


in Eq. ([5^ 


]) for multiplets 8, 10, 


and 10. 












J 


J(J + 1) 


C2(SUr(3))- J(J + l)-f 
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3/2 
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15/4 
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TABLE IL Expectation value of Eq. ([! 


MD, where S e 8, 10,T0. 
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8/12 


9/12 
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S G lO 
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6/12 




7/12 


8/12 


9/12 



where 

1 



^^ = i-J(J + ^)^ (22) 

3" 



.2 1 



C2(SUr(3))- J(J + i: ^ 



(23) 



K 

<1B = l{^f\l-D^:^{A)\^^^). (24) 

The quantities u and k distinguish the multiplets (8, 10, 10), and the values of J{J + 1) 
and C2(SUr(3)) — J( J + 1) — | are given at Table [B For the states with J = I, uj and R 
are regarded as the angular frequencies of the rotation in ordinary space and strangeness 
direction, respectively. The expectation value qb is a source of the SU(3) symmetry breaking 
on the profile function and characterizes each baryon state. Table HTl shows the values of qs 
for the individual 8, 10, and 10 baryons. Therefore, the profile functions derived from the 
classical Hamiltonian fl2T]) change the shapes according to the baryon states. Our approach 
modifies RRA in this meaning. 

We obtain the equation of motion for F{p) from Eqs. f|T9l) and fl2T]) : 

Cp.{p, F)F" + C(^,).(p, F) {F'f + CAP, F)F' + C(p, F) = 0, (25) 
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where 



Cf'{p,F) 
C{p, F) 



1 + -u^ sin^ F + -k^ sin^ - 

V3 4 2 
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3 4 2 



2 



+ -e'Co' ( 1 + 
4 



P / P 
2 sin^ F 



m^g sin F 



sin2F 

2 F 



sm 



l + (l + 3cosF) , 
p^ 



sinF, 



and rries is an effective meson mass given by 



(26) 
(27) 
(28) 



(29) 



(30) 



Coupled equations f|T2l) . f|T3|) . fl22l) . fl23|) . and fl25l) are self-consistently solved under the 
boundary conditions ([6]). Then, independent parameters are (e, mgfj) only. Therefore the 
effect of the flavor symmetry breaking on the profile function is expressed by rhes- In RRA, 
rfiefT = rries/^efT,) = m^/(e/^) ^ 1. In our approach, the effective mass can take a value 



niK according to qs- 



and niK = 495 MeV 



'his fact is important. If we estimate that e = 3.87, /tt = 44.5 MeV, 
31| . a large effective mass fhes ~ 2.8 is obtained. Therefore, we can 



expect qualitatively different behavior of the soliton solution in our approach. 



III. SOLITON SOLUTION IN THE MEAN FIELD APPROACH 

Our next task is to perform the self-consistent procedure in Sec. Ill B[ The procedure is 
faced with two kinds of complexity: the instability of the soliton solution and the dependence 
of the soliton solution on the multiplets. The instability results from Eqs. f l2S]) and (Q. The 
dependence on the multiplet is brought into the calculation by Eqs. and Therefore, 
we discuss these problems separately in the following sections. 
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A. Instability of the soliton solution 



To investigate the instability of the soliton solution, we treat (e^cj^, e^K^, m^fj) as input 
parameters in this subsection and the next one. The parameter space is designated as Ai. 

The stable soliton solutions of Eq. (1251) are obtained only in a restricted area of A4. The 
restriction has two origins. One is a behavior of the profile function F at p ~ oo due to the 



centrifugal force of the rotation 



F ~ (l + pp)e-^'', (31) 



p = ^lml,-e^(^l^^ + h^y (32) 

Here the rotation pushes F out of the center of the soliton. For the stable soliton solution, 
the following condition should be satisfied; 

ml, - + h^'j > 0. (33) 

Therefore, the rotating SU(3) Skyrmion with rhcs = is unstable. That is analogous to the 



result of the rotating SU(2) Skyrmion in the chiral limit [27|]. The rotating SU(3) Skyrmion, 
however, can exist in a limit (m^r = 0, rfiK 0), because m^^ > from Eq. flSCT]) and Table HTl 
Another origin of the restriction on is a behavior of the coefficient function of F" 
in Eq. fl2^ . We have the second condition for the stable soliton solution satisfying the 
boundary conditions ([6]): 

Cf.{p,F{p))>0. (34) 
For verification of this condition, we define a curve F* in (p, F) plane by 

Cf4p,F*) = 0. (35) 

Function F* should not be confused with the profile function F. We will show that F 
satisfying the boundary conditions ([HD cannot cross F* and this requirement is equal to 
condition flMD. 



At first, we investigate properties of F*. Since Eq. fl35|) is a quadratic equation for sin^ 
it has two formal solutions for fixed p and (e^a)^, e^k"^). These formal solutions, however, do 
not always support two real number values of F* in the range < F* < tt which has one-to- 
one correspondence with the range < sin^ < 1- If -F* is a real number solution in the 
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P 

FIG. 1: Typical forms of F* within -vr < _F* < vr. The cases of e^h^ < 4, e^K^ = 4, and e^K^ > 4 
are represented by the soUd Unes, the dashed hnes, and the dash-dotted hnes, respectively. For 
classification of the lines, we set e'^oj^ ~ 1. 

range, ±F* + 2mT {n: integer) also are solutions in other range. Practically, we can restrict 
the value of F* to the range — vr < F* < it, because the boundary conditions (jS]) ensure 
that the value of F is in the range. Figure [1] shows the typical forms of F* for e^K^ < 4, 
e^k"^ = 4, and e^k"^ > 4. The form of F* changes drastically at e^k'^ = 4. In particular, there 
are the constant solutions F* = ivr for e^K^ = 4, and F* for e^K^ > 4 always reaches p = 0. 
Figure [2] show also the dependence of F* on e^u)^ in the range < F* < vr. For a larger 
value of e^(D^, F* becomes closer to the axes p = and F* = 0. 

Next, we explain how to verify Eq. (!34|) . If F crosses F* at some radius p = p*, Eq. (!25|) 
becomes a quadratic equation for F\p*): 

C^p,^2{p\F*) {F'f + Cf'{p\F*)F' + C{p\F*) = Q. (36) 

Then we can statically calculate the value of the discriminant 

D* = CF'{p\F*f -AC^j„)2{p\F*)C{p\F*). (37) 

Of course, D* < means that Eq. fl36|) has no real number solution and F cannot cross F* 
at p = p* from the beginning. Moreover Eq. (125|1 does not have the real number solution in 
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e^K^ < 4 




fi4i,2 = 2 e^Cu^ = 1 



_±_ 



2 4 

e4;i2 = 4 



T" 



T" 



e^cD" = 3 




FIG. 2: Typical dependence of F* on e^u)^ in the range < F* < TT for e^K^ < 4, e'^R^ = 4, and 
e'^K^ > 4. Actual values of e'^K? in the figures are 3, 4, 5, respectively. 

a neighborhood of the point, because its discriminant for F': 



CAp.Ff - AC^^,^.{p,F)C{p,F) 

- 4C^^,f{p,F)CMp,F)F" 



(38) 



approaches the value of D* (< 0) near F*, and the value of F'{p) becomes complex numbers. 

For D* > 0, the values of F'{p*) are real numbers. However it is analytically unclear 
whether the values are consistent with the boundary conditions (Q. From numerical calcu- 
lations, we conclude also here that F and F* cannot cross. 
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P 

FIG. 3: Profile functions Fi {i = 0,1,2,3) and curves F* for e^k^ < 4 {e^K^ = 3, e^a)^ = 2, 
fhcs = 2). The solid lines represent F^. The dashed lines indicate D* < parts of -F*, the dash- 
dotted lines D* > parts, and marks "•" represent the points at which D* = 0. Marks "x" 
represent breakdowns of Eq. (j25p at these points. 

Figures [3] and m show the examples. Profile functions Fi {i = 0, 1, 2, 3) and F* for e'^k'^ < 4 
are plotted in Fig. [31 Each of Fi corresponds to different values of F'{0), namely 

> F^(0) > Fm > Fm > F[{0). (39) 

Only Fq satisfies the boundary conditions In the cases of Fi 2, there are breakdowns of 
Eq. fl25l) . because the profile function approaches the D* < part of F*. Profile function 
F3 survives but does not satisfy the boundary conditions due to the scattering by the 
D* > part. 

Figure m shows Fi {i = 0,1,2,3) and F* for e^/t^ > 4. Also here, only Fq satisfies the 
boundary conditions ([6]), and F12 break because of the same reason as that of the case 
e^K^ < 4. Moreover, there is a new situation that F* divides F3 from the other profiles (e.g. 
F2). Therefore, F2 is in Cp" > area and F3 is in Cp" < area. Although the values of 
^2(0) and -^3(0) are close, the profile functions F2 and F3 separate with an increase in p. It 
seems that there is a repulsion between F and the D* > part of F*. Profile function F3 
survives near the D* > part at small p but breaks near the D* < part at last. 

From these results, we conclude that F satisfying the boundary conditions cannot 
cross F*. Therefore, F* divides (p, F) plane into two areas: Cp" > area and Cp" < area. 
Profile function F{p) lives in only one area containing the boundary point (p, F) = (00, 0) 
of Eq. (|6]). Since Cpn {00, 0) = 1 > at the boundary point, we obtain condition (134]) of the 
stable soliton solution. 
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FIG. 4: Profile functions Fi {i = 0,1,2,3) and curves F* for e^K^ > 4 (e^K^ = 10, e'^cj^ = 0.5, 
rhes = 1.9). Meanings of the figure symbols are the same as those of Fig. [3l 

We should choose the value of F'{0) carefully so that F is away from the D* < part 
of F*. Since F* approaches the axes p = and F* = at the larger values of e^u)^ and 
e^K^ (Figs. [Hand [2]), the choice of F'{0) becomes more difficult. This situation improves for 
larger value of rrigfj, because F dumps faster according to Eq. fl3Tl) . Thus, with an increase 
in fhes, the area given by Eq. (IMll enlarges with rhes in parameter space Ai. 

B. Profile function of the stable soliton solution 

In this subsection, we discuss several constraints on the parameter space Ai and the 
stable sohton solutions of Eq. ( l25i) . At this stage, we have two conditions ( l33l) and ( IMl) for 
the stable soliton solutions. Condition (!33|) is explicitly parametrized by (e'^u;^, e'^/t^, rrigg), 
and condition flM|) implicitly. Both the conditions define together an area of the stable 
soliton in Ai. In Fig. [5], we show a critical surface that separates the areas of the stable 
solution and the unstable solution in A4. The area of the stable soliton is on the upper side 
of the surface. 

Condition ( l33i) defines the critical surface for small value of fhes, and condition does 
so for large value of meff. There is a boundary curve, at which these conditions change 
the roles on the critical surface. The boundary has rhes ~ 1-6. We designate the lower 
(higher) critical surface as Sl{h)- Sl is plane, because Eq. fl33l) is a linear equation for 
(e'^u)^, e^K^, rrigjj). Sh curves upward steeply. 

Here, we illustrate the effects of the rotation and the symmetry breaking with the four 
points R, A, B, and C depicted in Fig. [5l Point R is placed at (e^u)^, e^K^, rh-gg) = (0,0,4) 



13 



FIG. 5: Critical surface given by Eqs. (|25p and Q in parameter space M. The sizes of the mesh 
on the surface are (Ae^u)^, Ae^K^) = (0.5, 2.0). The area of the stable soliton is on the upper side 
of the surface. The solid line on the critical surface indicates the boundary curve (meff ~ 1.6) at 
which two conditions (j33p and (j34p change the roles. Point R is placed at (0,0,4). Points A, B, 
and C are placed at (3.5,0,4), (0,12.7,4), and (1.5,9,4) respectively, and are slightly above the 
surface. 

in A4. Points A, B, and C are placed at (3.5,0,4), (0,12.7,4), and (1.5,9,4) respectively, 
and they are slightly above the critical surface Sh- Figure [6] shows the profile function, its 
derivatives, and curve F* [Eq. fl5^ ] with the parameters corresponding to points R, A, B, 
and C in . 

Case R corresponds to RRA, because the influence of the rotation is ignored. In Fig. [6l 
there is a flat part of F" at p ~ 0.5. It is caused by the large meson mass [rhcs ~ 2). 
Although such a large meson mass is unfamiliar in other studies on the proflle function, it is 
legitimate in our approach as noted in Sec. Ill B[ If rhes increases more, the flat part dents 
downward. However, F" does not cross the zero, and F' monotonically increases. 

The proflle function with small rhes (the parameters near Sl) is given by spreading the 
solution of case R according to Eq. fl3Tl) . Then, the effects of the rotation and the symmetry 
breaking appear in the single mass parameter /i. If meg decreases further, /i becomes a 
complex number and the soliton becomes unstable because of the emission of the real meson. 

In the cases of A, B, and C, there are further characteristic behaviors of the proflle 
functions. Since points A, B, and C are close to Sh in A^, the corresponding proflle functions 
emphasize features of the rotation. 

Gradients |-F'(0)| in the cases of A, B, and C are small compared with that in the case 
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of R, because the rotation pushes the profile function out of the center of the sohton. Then 
F" should change largely for F to maintain the asymptotic form fl3T]) with the large mefr- 
Therefore, F" crosses the zero, and the behavior of F' becomes complex. 

Case A represents the profile functions deformed by the rapid spatial rotation {e^k^ < 4 
and large e^u)"^). Curvature F" change intensely at intermediate p region. This profile 
function is apparent at this stage, but it is excluded by the self-consistent procedure as 
noted in the next subsection. Moreover, such a profile function is physically unimportant; 
point A in Fig. [5] is on e^k^ = plane corresponding to the SU(2) Skyrmion and its effective 
mass is large {rhes ~ 2), however, physically rhcs = rriT^/^ef-,,) <^ 1 in this sector. 

Case B represents the profile functions affected by the rapid flavor rotation (e^K^ > 4 and 
small e'^u^). Curvature F" is already negative at p ~ 0. Since F* reaches p = 0, Eq. fl51]) 
reduces to 

4 ~ 2 

1 + 2F'{0Y - ^ > 0, (40) 

according to the boundary conditions ([6]). If rfigg decreases, F'(0)^ becomes too small and 
Eq. (HOD fails. 

Case C represents the profile functions affected by both the spatial and the flavor rotation 
(e^K^ > 4 and medium value of e^uj'^). Although F* reaches p = and Eq. PUI) is valid 
here too, the spatial rotation affects the profile function at p ~ 1 and F" changes greatly 
there. This spatial rotation has smaller angular frequency than that of case A. Therefore, 
the flavor rotation enhances the effect of the spatial rotation. 

C. Self-consistent soliton solution and the dependence of the classical soliton on 
the multiplets 

We are ready to study the self-consistent solution of the coupled equations Eqs. (IT^ . (IT^ . 
fl22]) . (1251) . ([25D, and the boundary conditions (jH]). The independent parameters reduce from 
(e^u)^, e'^K^, rrigg) to (e,meff), and the quantities (cD^, k^) are self-consistently determined for 
each muhiplet 8, 10, and 10. From Tableland Eqs. (|22D and (|23|), one can estimate that 

u\^^5xu\, (41a) 
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for 8 and 10, and 



~ 2 I ~ 2 



8 ' 



(42a) 
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(42b) 



for 8 and 10, and 




10 ' 



(43a) 



10 



(43b) 



for 10 and 10. 

Since the number of the independent parameters is two, the self-consistent solutions form 
surfaces for each multiplet in parameter space Ai. The surfaces are limited by the critical 
surface in Fig. [51 We call the surfaces "self-consistent surfaces" and show these in Fig. [3 
In addition. Fig. [8] shows the intersection lines between the critical surface and the self- 
consistent surfaces on the contour map of the critical surface. The self-consistent surface 
for 10 is away from the others because of Eqs. (HTil . Further the surface for 10 is closer to 
e'^u^ = plane than the surface for 8 because of Eqs. (jUj). 

These figures are useful for relating the profile functions with multiplets 8, 10, and 10. 
For example. Fig. ([7]) shows that the profile functions in case A in Fig. [6] are excluded, 
because point A is on e^k"^ = plane in and any self-consistent surfaces do not pass 
through this plane except for the m^g axis. However the profile functions in the cases of R, 
B, and C can be self-consistent solutions. 

Figure [8] shows that if the self-consistent solution is evaluated near Sh, always e^K^ > 4. 
Then the profile functions of 8 and 10 are similar to that of case B in Fig. [6l and the profile 
functions for 10 are similar to that of case C. 

Also in parameter space (e,meff), there are curves that separate the area of the stable 
soliton solutions and the area of the unstable ones in each multiplet. We call the curves 
"critical curves" and show these in Fig. [9l Every curve has an upward ledge at fhes ~ 1.6. 
The left area of the ledge is restricted by Eq. ( l33i) . and the right one by Eq. ( IMI) . These 
curves correspond to the intersection lines between the critical surface and the self-consistent 
surfaces in Fig. [71 One can use this figure to decide whether the adjustable parameters 
(e, /tt, m^r, mx) admit the stable soliton solution through Eq. f[30|) . 
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FIG. 7: Self-consistent surfaces for the multiplets 8, 10, and 10 in parameter space M. The 
dashed hnes represent the self-consistent surfaces. Three panels show these surfaces one by one 
in the order of 10, 8, and 10. The solid lines denote the critical surface in Fig. [5l The dashed 
lines on the critical surface represent the intersection lines between these two kinds of the surface. 
The self-consistent surfaces under the critical surface are spurious. They have been shown due to 
a limit of the ability of our graphic software. 



The critical curve should reach {e,fhes) = (0,0) because of Eq. ( l33l) . However, it is 
difficult to decide whether the solution of Eq. (1251) is stable for rhes ~ 0, because the 
instability for Th^s ~ appears at the large radius p ^ 1. Therefore, we show the curves 
only for rhes > 0.2. 

For rficR > 1.6, the areas of the stable soliton become narrow, because the moments of 
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FIG. 8: Intersection lines between the critical surface and the self-consistent surfaces for the 
multiplets 8, 10, and 10. The dashed lines represent the intersection lines. The dash-dotted lines 
represent the contour lines of the critical surface with step Am^g = 2. The solid line indicates the 
boundary curve in Fig. \5\ 

10 
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FIG. 9: Critical curves for multiplets 8, 10, and 10 in parameter space (e,meff). These curves 
separate the areas in which the soliton is stable or unstable. The horizontal axis rhes is different 
from meff/-F7r = mes/C^fw) of Ref. [2g|. 



inertia are the decreasing functions of rhcs in our approach, and as a result the 

stability condition flMl) with Eqs. fl22l) and fl23|) becomes severe for parameter e. 

From Eqs. fHTl) . and it is reasonable that the critical curve for 8 is higher 

than those for 10 and 10 through all value of rhcs- On the other hand, from Eqs. ([33D, flMl) . 
and (j43l) . it is not clear which critical curve for 10 and 10 is higher. Indeed, the curves for 
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10 and 10 change their relative height according to rhcs- In particular, the critical curve 
for 10 is lower than that for 10 in the area rhes > 1.6. It is explained as follows. The 
self-consistent solutions corresponding to the critical curves are obtained near Sh in Fig. [71 
and the solutions have e^K^ > 4. Therefore, condition (!34l) reaches p = 0. For a large 
effective mass, condition fl34|) is effective at p ~ because of Eq. fl3T|) and it reduces to 
Eq. fl4Up containing only k"^. Since k"^ for 10 is about three times larger than that for the 
other multiplets in Eqs. f l42b[) and f l43b[) . the area of the stable solitons for 10 is narrower 
than those for 8 and 10. 

The terms proportional to in Eqs. and (140|) originate from the term 
in Eq. fl2Tl) . This term y/^^K^ is regarded as the coupling between the rotation into the 
strangeness direction and the amplitude sin(F/2) in The amplitude is the zero mode 



fluctuation around the hedgehog profile in the flavor symmetry limit [25j, and it represents 
an intrinsic motion on the soliton. If this term becomes large compared with Mq in the mean 
field Hamiltonian fl2Tl) . the collective rotation and the intrinsic motion cannot dynamically 
separate and RRA fails. That is a limit of the applicability of RRA pointed out from a 
general argument in Ref. 17|]. In our approach, the influence of this coupling is dynamically 
included in the calculation of the profile function through the mean field Hamiltonian ( l2Ti) . 
Therefore, condition represents a more realistic limit so that the coupling does not 
destroy the soliton itself. 



IV. BARYON MASS 

With the profile function of the stable soliton, the classical soliton mass (fTTI) . the moments 
of inertia f[T2]) . f[T^ . and the symmetry breaking (HM are evaluated. Then the eigenvalue 
equation fll6l) is numerically solved by the Yabu and Ando method, and the baryon masses 
are obtained. We introduce 

= / J2 (^Ebn - AE^^^f, (44) 
y Be8,io 

where AEbn is a difference in the predictive value of the mass between the baryon B and 
the nucleon N, and AE"^^ is its experimental value. An isospin multiplet is represented as 
a baryon B in this formula, because baryons in an isospin multiplet are described by the 
same soliton solution. The quantity AE'^" measures an error of the predicted baryon mass 
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splitting for the multiplets 8 and 10. 

The parameters (e, rhes) of the stable solitons are placed below the critical curves in 
Fig. O A parameter set {e, 771^^,771 k) corresponds to the 12 points in the parameter space 
(e, TTLcs) according to Eq. fl5U]) and Table [ITl While the Skyrme parameter e determines 
the vertical positions of the points, the masses give the horizontal positions and the 

spreads of the points. We discuss only parameters (e, rfiTr, mi^) which admit the existence 
of the 8, 10, and 10 baryons. Since {e,mT,,mK) are dimensionless, the energy scale {fn/e) 
cannot be specified, and the stability of the solitons is not sufficient to determine the values 
of all parameters (e, /^r, m^r, mi^-). Using the degrees of freedom, we fit the — A mass 
difference or the absolute value of the mass to its experimental value. 

Equation f[T7|) gives the experimental value of the N — A mass difference and the accurate 
baryon mass splitting (small AE'^"). The parameter set and the predicted baryon masses 
are shown as (1) in Tables UTTl and HV] respectively. Here, the value of /^r, 93 (MeV), is 
given by hand. A larger value of /^r gives a slightly smaller value of AE^^^, but it leads to 
extremely larger baryon masses in proportion to /^r; for example, AE^" = 74 (MeV) and 
En = 5992 (MeV) for e = 2.8, U = 186 (MeV), = 106 (MeV), and tuk = 747 (MeV). 
Table IVl shows the dimensionless effective mass {rhes), the classical mass (Mq), the moments 
of inertia (a^, (3"^), and the symmetry breaking (7) for parameter set (1). In our approach, 
these quantities vary according to the baryon states: rhes characterizes each baryon state, 
Mq is the increasing function of rfiefr, and (a^, jS'^, 7) are the decreasing functions. 

However, Eq. ( |T7I) cannot give the experimental value of the A^ mass and the accurate 
mass splitting simultaneously in our approach. In particular, the values of the baryon masses 
are lar ge c ompared with the observed ones. That is a common phenomena in the Skyrme 



model 



3l| . but the tendency is more severe in our approach. The experimental values of 



the baryon masses lead to the smaller value of which is the energy scale in this model. 
Therefore the value of the Skyrme parameter should be large for the rotational energy to 
generate the mass splitting. However, since the Skyrme parameter is restricted by the critical 
curves for 10 and 10 in Fig. [9], the magnitude of the mass splitting is not sufficiently large. 



There is a more fundamental method 
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36|, 



371, 



38l | for the Skyrme model to reproduce 



the observed baryon masses. In this method, the Casimir energy (< 0) due to the existence 
of the soliton is added to the baryon masses. In RRA, the Casimir energy is the quantity of 
order A"^, and the value is common to all baryon states. Therefore the addition of this energy 
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does not change the mass sphtting, and one can discuss the mass sphtting and the values of 
the masses separately. In our approach, the Casimir energy changes its value according to 
the baryon states as well as the shape of the soliton and contributes to the mass splitting too. 
Thus, our self-consistent procedure should include the effect of the Casimir energy to treat 
the mass splitting and the masses themselves simultaneously. However that is a complicated 
task to be examined in detail here, because the simple analytic form of the Casimir energy 
is not known. 

Instead, we adopt a subtraction method 3l|] to estimate the effect. In this method, the 
unsubtracted mass formula (ITTI) is replaced by the subtracted one: 

% + ^S^SB~S,l (45) 



3/?2 2/32 

where the quantity £q is the lowest eigenvalue of Eq. (fTSj) corresponding to the vacuum- 
like state with (/, J, "K, Yr) = (0,0,0,0). Equation fH51) improves the behavior of Eb by 
removing the vacuum fluctuation energy according to the increase of the symmetry breaking 
and reproduce the mass splitting accurately [31 1. 

Table IIVI shows the baryon masses calculated by Eq. fH5l) with parameter sets (2) and 
(3) given at Table IIIII Set (2) fits the N — A mass difference, and set (3) fits the mass. 
Both the parameter sets give the accurate mass splitting. In addition. Table I VI I gives the 
values of meg, Mq, a^, P'^, and 7 for set (2), and Table rvTTl gives those for set (3). 

The deformation of the soliton reproduces the mass splitting accurately for any param- 
eter set given at Table IIIII and it has the sizable effects on the 8, 10, and 10 baryons 
masses as seen from Tables |Vl IVIt and IVIII However the mass splitting is caused by the 
different terms of the Hamiltonian according to the parameters. For example, the con- 
tributions of these terms to the N — E mass difference are estimated at the difference 
in (Mq, ^a^w^, 1/3^^2, It^b) calculated to S and A^. In RRA the N — E mass differ- 
ence is dominated by the symmetry breaking term. However, in our approach, it is dis- 
tributed as follows: (147,7,60,126) (MeV) for set (1), (51,11,72,254) (MeV) for set (2), 
and (17,14,98,201) (MeV) for set (3). Therefore the effects of the rotation and the sym- 
metry breaking mix through the deformation of the soliton each other, and the sizes of the 
effects are large. 
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TABLE III: Parameter sets. "Exp." denotes the experimental values. 



Set 



(2) 



(3) 



Exp. 



fn (MeV) 
(MeV) 
rriK (MeV) 



3.00 
93 
196 

1042 



3.40 
147 

616 



6.05 

46 
56 
551 



93 
140 
496 



If the baryon masses are given by parameter set (3), the 10 and 10 baryons are affected 
obviously by the critical curve, because the set corresponds to the points (e = 6.17, rhcs = 
1.36 — 1.81) just below the critical curves for 10 and 10 in Fig. O Then, as mentioned in 
Sec. nil CI the profile functions of 8 and 10 are similar to that of case B in Fig. [6], and the 
profile functions of 10 are similar to that of case C. 



V. SUMMARY 

In this paper, we have investigated the profile function of the SU(3) Skyrmion depending 
on the octet, decuplet and antidecuplet baryon states. The equations of motion for the 
profile function are given by the variation of the mean field Hamiltonian. The Hamiltonian 
is the expectation value of the collective Hamiltonian operator for the baryon state and 
depends on the profile function itself through the moments of inertia. Thus, we should solve 
the equations of motion self-consistently. As the result, the profile function is affected by 
not only the rotation of the Skyrmion but also the flavor symmetry breaking. 

The influence of the symmetry breaking on the proflle function is represented by an 
effective meson mass which varies according to the baryon states. The effective mass in the 
rigid rotator approach is the pion mass and usually small. In our approach, the effective 
mass can take a value about the kaon mass and, the qualitatively different behavior of the 
soliton solution emerges. 

In general, the rotation pushes the proflle function out of the center of the soliton, and 
the symmetry breaking (the effective meson mass) attracts the proflle function. For small 
effective mass, these effects are represented by the single mass scale in the asymptotic form of 
the proflle function. Then the instability of the soliton appears as the spontaneous emission 
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TABLE IV: Baryon mass differences from the nucleon mass for parameter sets (1), (2), and (3). 
"Exp." denotes the experimental values. AE^'^'^ is defined by Eq. (j44p . Only row of N gives the 
absolute values of the nucleon mass. Marks "*" denote the input values for the energy scales. All 
units are (MeV). 

Set (1) (2) (3) Exp. 

A^;^'''' 90 104 119 

8 

iV (abs.) 3618 3483 939* 939 

A 187 177 154 183 

S 314 310 270 256 

E 396 382 326 379 
10 

A 293* 293* 307 293 

S* 445 435 442 445 

E* 570 554 551 595 

n 672 654 641 733 
To 

e+ 708 483 441 601? 

888 661 613 ? 

Si^ 1022 800 746 ? 



10 



1059 865 819 ? 



of the real meson and restricts the parameter space of the self-consistent solutions. That is 
similar to the result of the SU(2) Skyrmion. 

For large effective mass, the influence of the rotation appears in the three cases. First, 
the rapid spatial rotation leads to the large variation of the curvature {F") of the profile 
function at the intermediate radius. Secondly, the rapid flavor rotation leads to the negative 
value of the curvature at the small radius. Thirdly, the flavor rotation enhances the effect 
of the spatial rotation and leads to the large variation of the curvature at the intermediate 
radius. Although the first case is excluded by the self-consistent calculation, the last two 
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TABLE V: Dimensionless effective mass ttIcs and parts of the Hamiltonian for the parameter 
set (1): classical mass Mq, moments of inertia (a^, and symmetry breaking 7. 







(MeV) 
Mo 


a2 


(10~^/MeV) 


(MeV) 

/ 


8 












N 


2.60 


2707 


6.85 


2.10 


1628 


A 


2.92 


2770 


6.48 


1.95 


1393 




3.22 


2828 


6.19 


1.84 


1226 




3.35 


2855 


6.07 


1.79 


1159 


10 












A 


2.89 


2742 


6.89 


2.08 


1531 


S* 


3.07 


2779 


6.69 


2.00 


1406 


"* 


3.25 


2813 


6.51 


1.93 


1302 


n 


3.42 


2845 


6.36 


1.87 


1216 


To 












e+ 


2.69 


2647 


7.90 


2.45 


2087 




2.89 


2681 


7.63 


2.34 


1890 




3.07 


2714 


7.42 


2.25 


1732 


"* 

"To 


3.25 


2746 


7.23 


2.17 


1602 



cases can be the self-consistent solutions and restrict the parameter space. 

The independent parameters of the self-consistent solution are the Skyrme parameter and 
the effective meson mass. There are areas of the independent parameters allowed for each 
multiplet. The allowed value of the Skyrme parameter for the octet baryons is the largest for 
all value of the effective mass, and those for the decuplet and antidecuplet baryons change 
the relative size according to the effective mass. At the large effective mass, the allowed 
value of the Skyrme parameter for the antidecuplet baryons is smaller than those for the 
octet and decuplet baryons. 

The baryon masses are evaluated by the unsubtracted mass formula and the subtracted 
one, respectively. Then, the deformation of the soliton reproduces the baryon mass splitting 
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TABLE VI: Dimensionless effective mass and parts of tlie Hamiltonian for the parameter set (2). 





I'll <r 
'"(■11 


(MeV) 


2 

n 


flO~^/MeV) 


(MeV) 














N 


0.84 


3200 


5.71 


2.08 


1134 


A 


0.95 


3221 


5.33 


1.91 


990 




1.05 


3241 


5.04 


1.78 


886 




1.10 


3251 


4.92 


1.73 


843 


10 












A 


0.94 


3199 


5.88 


2.15 


1181 


S* 


1 on 


3211 


5 fiB 


2.04 




"* 


1 07 


3222 


5 4fi 


1 9fi 


1018 


Q 


1.12 


3234 


5.30 


1.89 


955 


To 












G+ 


0.87 


3174 


7.23 


2.74 


1693 




0.94 


3180 


6.88 


2.58 


1550 




1.00 


3187 


6.59 


2.46 


1434 


"* 
"To 


1.07 


3194 


6.35 


2.35 


1338 



accurately with both the mass formulas and has the sizable effects on the baryon masses. 
Therefore the effects of the rotation and the symmetry breaking cannot separate clearly. 

The subtracted mass formula can reproduce not only the mass splitting but also the 
observed masses, though the pion decay constant is too small. Since the formula is inspired 
by the Casimir effect, the Casimir energy should be investigated for our self-consistent 
procedure to solve the problem of the small pion decay constant. It remains as a matter to 
be researched further. 

Other physical properties (e.g. magnetic moment, charge radius, etc. . . ) are affected by 
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TABLE VII: Dimensionless effective mass and parts of the Hamiltonian for the parameter set (3). 





''^eii 


(MeV) 


o? 


(W~^ /MeV) 

0^ 


(MeV) 














N 


1.36 


578 


4.77 


1.85 


917 


A 


1.54 


584 


4.4 


1.67 


776 




1.70 


591 


4.14 


1.54 


675 




1.77 


596 


4.05 


1.49 


635 


10 












A 


1.52 


606 


5 31 


2 05 


1 01 3 


y;* 




621 


5 1 


1 94 


90? 


"* 


1 72 


640 


4.97 


1.84 


808 




1.81 


660 


4.88 


1.76 


731 


To 












G+ 


1.41 


613 


6.87 


2.84 


1586 




1.52 


619 


6.49 


2.66 


1433 




1.62 


627 


6.20 


2.52 


1307 




1.72 


640 


5.98 


2.41 


1198 


the deformat 


ion. The 


study in this direction 


is in prog 
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